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We propose a theoretical framework for predicting the protocol dependence of the jamming tran¬ 
sition for frictionless spherical particles that interact via repulsive contact forces. We study isostatic 
jammed disk packings obtained via two protocols: isotropic compression and simple shear. We show 
that for frictionless systems, all jammed packings can be obtained via either protocol. However, the 
probability to obtain a particular jammed packing depends on the packing-generation protocol. We 
predict the average shear strain required to jam initially unjammed isotropically compressed pack¬ 
ings from the density of jammed packings, shape of their basins of attraction, and path traversed 
in configuration space. We compare our predictions to simulations of shear strain-induced jamming 
and find quantitative agreement. We also show that the packing fraction range, over which shear 
strain-induced jamming occurs, tends to zero in the large system limit for frictionless packings with 
overdamped dynamics. 

PACS numbers: 45.70.-n,61.43.-j,64.70.ps,83.80.Fg 


I. INTRODUCTION 

Dry granular materials are composed of macro-sized 
particles that interact via repulsive contact forces. Due 
to dissipative interactions between grains, granular ma¬ 
terials exist as static packings in the absence of exter¬ 
nal driving [l|. As a consequence, granular packings are 
out-of-equilibrium, and their structural and mechanical 
properties depend on the protocol used to generate them. 
Experimental packiM-generation protocols includ e g ravi- 
tational deposition , vibration [3| , compression , and 
shear Several computational studies have pointed 

out that the distribution of jammed packing fractions de¬ 
pends on the compression rate [3, Q and rate at which 
kinetic energy is removed from the system I, [i3. In 
addition, experimental studies of photoelastic disks have 
identified key differences between granular packings gen¬ 
erated via isotropic compression and pure shear |ll| . 

There has been significant work on understanding the 
scaling behavior of the elastic moduli and contact num¬ 
ber near jamming onset in model granular packings com¬ 
posed of frictionless spherical particles generated using 
isotropic compression [l2|. However, there is currently no 
theoretical understanding of how the ensemble of static 
packings and their properties vary with the protocol used 
to generate them. For example, what is the difference 
in the distribution of jammed packings generated via 
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isotropic compression versus shear? 

We focus on isostatic jammed packings of frictionless 
disks generated via different combinations of isotropic 
compression and simple shear and study the distribu¬ 
tion of jammed packings as a function of the path taken 
through conhguration space. A recent study has dis¬ 
tinguished between ‘compression-only’ jammed packings 
that possess non-zero pressure and positive shear moduli 
for some but not all boundary deformations, and ‘shear- 
stabilized’ jammed packings that pos sess positive shear 
moduli for all shear deformations [T^- We describe the 
protocol dependence of compression-only jammed pack¬ 
ings, which are experimentally realizable [14| and are rel¬ 
evant for understanding jamming in systems with fric¬ 
tional interactions [T5j |. 

We find several important results. First, an exponen¬ 
tially large but finite number of jammed packings with 
an isostatic number of contacts Nc = = 2N' — 1 

(where N' is the number of disks in the force-bearing 
backbone) exist in configuration space, defined by the 
disk positions, packing fraction (p, and shear strain de¬ 
formation 7 of the system boundaries. In small systems, 
nearly all jammed packings can be enumerated [l4l| . For 
example, we have shown that isostatic jammed packings 
form one-dimensional geometrical families as a function 
of shear strain [T^. We will show that the choice of the 
packing-generation protocol does not change the ensem¬ 
ble of isostatic, jammed packings, but instead changes 
which packings are visited during particular trajecto¬ 
ries through configuration space. The average properties 
of jammed packings change for different protocols be- 
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FIG. 1: Fraction of jammed packings (or cumulative distri¬ 
bution) C{(j>, 7 ) along different paths in the packing fraction (f) 
and shear strain 7 plane for N = 32. (a) 7 ) for protocol 

1 (isotropic compression at fixed 7 ; solid line) and C 2 (<(>, 7 ) for 
protocol 2 , i.e. compression to (f) followed by shear strain to 
7 = 0.1 (dotted line), 0.3 (dot-dashed line), and 0.5 (dashed 
line). I, II, and III indicate the packing fractions in (b). (b) 
We show (dashed lines) and C 2 (<)>, 7 ) (solid lines) at 

fixed (f> = 0.815, 0.824, and 0.832 indicated by I-III. Proto¬ 
col dependence can be seen in the difference between Ci and 
C 2 evaluated at the same (j> and 7 , e.g. a.t cj) = 0.824 and 
7 = 0.67 as highlighted by the dashed double arrow. Right 
and left solid arrows indicate protocols 1 and 2 , respectively. 


cause the probabilities for obtaining each jammed pack¬ 
ing varies with protocol. 

We develop a theoretical description of the protocol 
dependence of the distribution of jammed disk packings 
(Fig. [T]). We assume that an initially unjammed system 
will jam when it encounters the basin of attraction of a 
jammed packing as it travels through configuration space. 
The probability to obtain a jammed packing is deter¬ 
mined by two factors: 1) the number density of jammed 
packings in configuration space, which is independent of 
the packing protocol, and 2) the path traveled through 
configuration space, which depends on the protocol. Us¬ 
ing this framework, we predict the average shear strain to 
jam initially unjammed packings at each (j) and show that 
the predictions agree with simulations of shear strain- 



FIG. 2: (a) Packing fraction (j) versus shear strain 7 for all 
isostatic jammed N = 6 disk packings. The solid black line 
obeys (j> = 71(7 “ 7 o)^ + <(>o with A = 0.776, (po = 0.665, 
and 70 = 0.35. Filled circles (downward triangles) indicate 
packings with positive (negative) local slope. The solid verti¬ 
cal arrow indicates protocol 1 , and the dashed vertical arrow 
followed by the dashed horizontal arrow indicates protocol 2 
that was used to reach a jammed packing at 7 = 0.8 and 
(j) = 0.725. (b) Jammed packing fraction <(>( 7 ) using protocol 
2 at fixed if) in the range 0.64 < (f> < 0.77 for N = 6. (c) Num¬ 
ber of = 32 jammed packings at each (j> and 7 (increasing 
from dark to light) from protocol 1. 


induced jamming. Our results indicate that the packing 
fraction range, over which shear strain-indnced jamming 
occurs, vanishes in the large-system limit for overdamped 
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frictionless systems. 

The remainder of the manuscript is organized into 
three sections. In Sec. m we describe our simulation 
methods and introduce the two protocols used to gen¬ 
erate isostatic jammed packings. In Sec. uni we show 
results concerning the protocol dependence of the distri¬ 
bution of jammed packings. We then describe a theoret¬ 
ical model that allows us to calculate the probability to 
obtain jammed packings as a function of the path that 
the system traverses in configuration space. In Sec. IIVI 
we summarize our results and conclusions. 


II. METHODS 

We study systems containing N frictionless bidisperse 
disks in a parallelogram with height L = 1 in two dimen¬ 
sions (2D) that interact via purely repulsive linear spring 
forces with energy scale e [I 4 ] . (Studies of bidisperse fric¬ 
tionless spheres in three dimensions are included in Ap¬ 
pendix]^) The mixtures contain half large and half small 
particles with mass m = 1 for both and diameter ratio 
ctlIc^S = 1-4. We employ Lees-Edwards simple shear- 
periodic boundary conditions, where the top (bottom) 
images of the central cell are shifted to the right (left) by 
TyL and 7 is the shear strain [T^ . We varied the system 
size from = 6 to 512 particles. 

Below, we describe results for two protocols to gener¬ 
ate jammed packings in the ^-7 plane. (See Fig. [5] (a).) 
Protocol 1 involves isotropic compression at fixed bound¬ 
ary shape parametrized by shear strain 7 . We start with 
random initial disk positions at (j)Q < 0.5. We succes¬ 
sively compress the system by increasing particle radii 
uniformly in small packing fraction increments dcj) and 
minimize the total potential energy per particle V/{Ne) 
(at fixed 7 ) after each step. Jamming onset occurs when 
k^nax > V/{Ne) > 0, with t4iax = 10“^®, or an equiv¬ 
alent threshold on pressure. For protocol 2, we start 
by isotropically compressing systems (at 7 = 0 ) to 4>, 
and if the system is unjammed with V/{Ne) <C Idnax, 
we successively apply simple shear strain to each parti¬ 
cle x'^ = Xi + d"fyi in strain steps dy < 10 “® followed 
by minimization of V/{Ne). We then identify the to¬ 
tal shear strain 7 at which the system first jams with 
10“^® > V/{Ne) > 0. Protocols 1 and 2 generate 
compression-only jammed packings. In Appendix lUl we 
also describe results for a third protocol, which is simi¬ 
lar to protocol 2 , but with simple shear replaced by pure 
shear. Note that isostatic jammed packings can also be 
generated using stress-controlled packing-generation pro¬ 
tocols [llj. 

III. RESULTS 

We display the cumulative distributions Ci 7 (()), 7 ) of 
jammed packings from protocols 1 and 2 in Fig.jTJ In (a), 
we show that applying shear strain increases the fraction 


of jammed packings at each (j), i.e. 02 ( 7 , (j)) shifts to lower 
(j) with increasing 7 . In (b), we show Ci(^, 7 ) obtained 
via isotropic compression versus boundary shape for (j) = 
0.815, 0.824, and 0.832 (corresponding to Ci{cj),0) Ki 0.2, 
0.5, and 0.8). I — Ci(^, 0) of the packings are initially 
unjammed at 7 = 0 and (j), and these jam with increasing 
7 as shown by the solid lines. By combining different 
amounts of shear strain and isotropic compression, the 
fraction of jammed packings at a given (j) can be tuned 
over a wide range, e.g. from 0.2 to 0.8 for packings at ^ = 
0.815. These results emphasize that the distribution of 
jammed packings depends strongly on the path through 
configuration space, e.g. protocols 1 and 2 indicated by 
the arrows in Fig. [T] (b). 

To understand protocol dependence, we examine the 
distribution of jammed packings in the t/f-y plane. In 
Fig. UK a), we show the packing fraction at jamming on¬ 
set 4> versus 7 for = 6 from protocol 1 (solid vertical 
arrow in Fig. [2] (a)). We find several striking features. 
First, jammed packings occur as geometrical families 
{i.e. segments of parabolas that correspond to jammed 
packings with the same interparticle contact network) 
in the ^-y plane. For N = 6, we are able to enumer¬ 
ate nearly all geometrical families over the full range of 
7 [ 1 ^ . When straining an initially unjammed system 
toward positive 7 at fixed (j) {e.g. horizontal arrow in 
Fig. [5] (a)), it will jam on a geometrical family with neg¬ 
ative slope {—\d(j)/dj\). For negative slopes, continued 
shear strain leads to overcompression, whereas for posi¬ 
tive slopes, continued shear strain leads to unjamming. 
This behavior is shown in Fig.[2](b) for protocol 2 at fixed 
ij) in the range 0.64 < </) < 0.77 for N = Q. Note that any 
of the jammed packings in Fig.[5](b) from protocol 2 and 
defined by {K}, (f’^ 7 can be generated using protocol 

1 with initial condition {K} and boundary deformation 
7 . As a result, we can generate the same jammed pack¬ 
ing at a given (j) and 7 using different combinations of 
compression and shear strain. We find similar behavior 
to that in Fig. [2] (a) and (b) for larger N, except that 
the parabolic segments in (j){'y) become smaller and more 
numerous, and thus geometrical families more densely 
populate configuration space (Fig.[2](c)). The number of 
jammed packings (at a given (j) and 7 ) becomes indepen¬ 
dent of 7 for A^ > 32. 

We develop a theoretical model using an analogy with 
absorption problems to calculate the probability to ob¬ 
tain isostatic jammed packings as a function of the path 
that the system traverses in configuration space. In prin¬ 
ciple, the number density of jammed packings T depends 
on the 2N coordinates of the disks, (j), and 7 , but not the 
packing-generation protocol. After integrating over the 
2N disk coordinates, is a function of 4> and 7 . However, 
we assume that the number density T{4>) is only a func¬ 
tion of (/) since T becomes independent of 7 in the large- 
N limit (Fig. [3] (c)). We imagine that a one-dimensional 
trajectory /!(<(), 7 ) through configuration space will en¬ 
counter the basin of attraction |19| of a jammed pack¬ 
ing with a probability T{(j))S{(j))dC during a step of size 
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FIG. 3: (a) Natural logarithm of the fraction of un¬ 

jammed packings (normalized by the (^-dependent decay fac¬ 
tor, J^{<j))S 2 {(l>)£ 2 { 4 >)), during protocol 2 at fixed ij) in the range 
0.798 < (j) < 0.844 (solid lines) for N = 32. The dashed 
line has slope —1. (b) Comparison of —dln[M 2 (((), 7 )/Mo]/d 7 
(open circles) from protocol 2 and —dln[Mi{(l))/Mo]/d(j> from 
protocol 1 with S 2 { 4 >) oc (solid line) or S 2 i 4 >) oc 

Si{(j>) (dashed line) for N = 32. The inset shows the dis¬ 
tances iiitj)) and a£2{4>) (with a ~ 5.5) traversed in conhg- 
uration space during protocols 1 (solid line) and 2 (dashed 
line) for N = 32. 


dC in configuration space, where S{(j)) is the 2N — 1- 
dimensional cross-section of the basin of attraction of a 
jammed packing perpendicular to dC. 

Thus, for protocol 1 with trajectories only along the 
decrease in the number of unjammed packings dMi{(j)) 
(or equivalently, increase in the number of jammed pack¬ 
ings) during a compression step dcj) is 

dMi(</)) = -Mi{^)T{cl2)Si{^)£i{^)d^, (1) 

where dC = £i{<j))d<j), £i{4>) is the distance in configura¬ 
tion space traversed during step dcj) at cf), and Sijcf)) is 
the average cross section for protocol 1. Eq. [T] can be 
solved for the number of unjammed packings at cj) during 
protocol 1 : 


Mi{cj)) = Mq exp 


- / F{ci>')SM)h{cl,')dcl,' 


'00 


, ( 2 ) 


FIG. 4: (a) Average shear strain 7 j(<(>) to jam an initially 
unjammed configuration at 4> and 7 = 0 (protocol 2) for N = 
32 (circles), 128 (triangles), and 512 (squares). We compare 
7 j((/!>) from protocol 2 to 1/{T{cj))Si{cj))£\{cj))) (Eqs. [3] and |4]) 
from protocol 1 for the same system sizes: N = 32 (solid line), 
128 (dashed line), and 512 (dotted line), (b) Distribution of 
jammed packing fractions Vijcf)) from protocol 1 for N = 32 
(solid line), 128 (dashed line), and 512 (dotted line) compared 
to predictions obtained from Vijcj)) = —MQ^dMi{cj))/dcj) with 
Mi{(j)) given by Eq.[5]and T{(j>)Si{(j>)li{(j)) given by Eq. fusing 
the measured value of 7 ^ for N = 32 (circles), 128 (triangles), 
and 512 (squares). 


where Mq is the number of unjammed packings at (j)o. 
For protocol 2, with trajectories only along 7 , we obtain 
a similar expression for the number of unjammed pack¬ 
ings: dM 2 {(j),"i)/d'y = —M 2 ((j),"i)£F{cj))S 2 {(j))£ 2 {cj)), where 
^ 2 ( 0 , 7 ) = Mi{(j)) eyi-p[-T((j))S2(,cl))l2{4>)l\, S 2 i^) is the 
average cross section for protocol 2 , and £ 2 ( 4 ’) is the 
distance traversed in configuration space for each shear 
strain step dj. 

Fig. [3] (a) shows that the fraction of unjammed pack¬ 
ings M 2 {(j),j)/Mi{cj)) decays exponentially with 7 during 
protocol 2 at each cj). This result emphasizes that we can 
calculate the product J-{(j))S((j))i{(j)) without enumeration 
of all jammed packings by measuring the decrease in the 
number of jammed packings with increasing shear strain. 
In the zeroth order approximation, the cross section Sjcj)) 
is independent of the path in configuration space and the 
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distance i traveled during each dcj) or dj step is constant. 
In Fig.[3](b), we compare T{(j))S{(j))£ from protocol 2 with 
the similar quantity —d\n.[Mi{(j))/MQ]/d4> from protocol 
1 and find qualitative agreement. 

We then independently measured (</>), defined by 
the accumulated distance in configuration space between 
the initial packing at cj) and relaxed packing at (j) + dcj), 
for protocol 1. We performed similar measurements for 
^ 2 ( 0 ), which gives the accumulated distance in config¬ 
uration space between the initial packing at 7 and re¬ 
laxed packing at j + dj for protocol 2. We show that 
the two are proportional to each other, £i{(j)) = a£ 2 {(j)) 
with a ~ 5.5, in the inset of Fig. [3] (b). By calculat¬ 
ing —d\a[Mi{(j))/MQ\/d(j) = F{4>)Si{4>)£i{4>) for protocol 
1 , we can compare to the similar quan¬ 

tity, !F{(j))S2{4’)£2{4’)^ for protocol 2. In this case, we 
assume that the cross section depends on the path in con¬ 
figuration space, e.g isotropic compression increases the 
overlaps of all interparticle contacts, while shear strain 
increases some but decreases others. In Fig. [3] (b), we 
show excellent agreement for !F{(j))Si^ 2 {(l))£i, 2 { 4 ’) for pro¬ 
tocols 1 and 2 for N = 32 provided we assume that 
S' 2 ((/>) oc Si{(j))£i{(p), and find similar quantitative agree¬ 
ment for all system sizes studied. (Additional details of 
the theoretical model are included in Appendix [Dll In¬ 
dependent measurements of Si^ 2 i 4 >) be performed in 
future studies. 

We now use the theoretical description of the protocol- 
dependent probability to jam to predict the average shear 
strain required to jam an initially unjammed isotropically 
compressed configuration at (j) and 7 = 0 : 

- - r _ 1 

Jo ^ ^ ^ 

In Fig. m (a), we show that the prediction for 
obtained from measurements of £F{(j))Si{(j))£l{(j)) using 
isotropic compression, agrees with simulations of shear 
strain-induced jamming. We find that 7 ^- grows rapidly 
with increasing system size and only packings with (j) > 
0.84 are jammed in the large-system limit [^. (We find 
similar results for applied pure shear in Appendix [Cj) We 
also calculate the distribution 'Pi{(j)) of jammed packing 
fractions (for isotropic compression) using data from pro¬ 
tocol 2. In Fig. |4](b), we show that Vi{4>) from protocol 
1 and Vi{(j)) = —MQ^dMi{(j))/d(f) with Mi{4>) given by 
Eq.[2]and J^((/))S'i((/))£i(((>) given by Eq.|4] (using the mea¬ 
sured value of 7 j ) collapse for all system sizes studied. 
The width of for isotropic compression narrows 

as 1/7V^ with A Ri 0.55 ± 0.05 and the peak approaches 
0.84 in the large-system limit (^ . 

IV. CONCLUSION 

In this manuscript, we developed a theoretical descrip¬ 
tion for jamming onset that allows us to predict the 


fraction of isostatic jammed packings that occur at (j) 
and 7 in terms of the path traversed in configuration 
space. This framework provides predictions for the av¬ 
erage shear-strain required to jam initially unjammed 
packings produced by isotropic compression, which agree 
quantitatively with simulations of strain-induced jam¬ 
ming in two- and three-dimensional systems subjected to 
simple and pure shear. In particular, we showed that the 
packing fraction range, over which strain-induced jam¬ 
ming occurs, shrinks to zero in the large-system limit for 
frictionless systems with overdamped dynamics. In fu¬ 
ture studies, we will investigate the role of static friction 
in stabilizing strain-induced jamming of dilute granular 
packings 
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Appendix A: 3D bidisperse packings 

In this section, we present our studies of compres¬ 
sion and shear-strain induced jamming of 3D bidisperse 
spheres, which are qualitatively similar to the results for 
2D bidisperse systems presented in the main text. We 
presented results in the main text on the distribution of 
jammed packing fractions V{(j)) and average shear strain 
7 j required to induce jamming in an originally unjammed 
configuration for systems composed of bidisperse disks in 
two spatial dimensions {d = 2). However, these results 
apply more generally than simply to two-dimensional 
packings of disks. The theoretical analysis in the main 
text described trajectories in the dV-dimensional config¬ 
uration space in which jammed packings exist, where d is 
the spatial dimension. The 2N-dimensional configuration 
space is already large, and thus we expect qualitatively 
the same results for three-dimensional (3D) sphere pack¬ 
ings, which exist in a configuration space that is only 50% 
larger, as we found for two-dimensional systems. 

We studied systems containing N frictionless bidis¬ 
perse spheres in a parallelepiped with sides of length 
L = 1 that interact via purely repulsive linear spring 
forces. The bidisperse mixtures contain half large and 
half small particles, both with mass m = 1 , and diameter 
ratio (Tl/(Ts = 1.4. As in 2D, we employ Lees-Edwards 
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simple shear-periodic boundary conditions, where the top 
(bottom) images of the central cell are shifted to the right 
(left) by 7 L, where 7 is the simple shear strain. 



(t> 


FIG. 5: Distribution of jammed packing fractions Pi(0) 
from protocol 1 for = 64 bidisperse spheres (solid 
line), compared to predictions obtained from = 

— MQ^dMi{(j))/d(l) with Mi(0) given by Eq. [2] and 
T{(f>)Si{(f>)h{4>) given by Eq. fusing the measured value of 
7 j for A'’ = 64 bidisperse spheres (circles). 


Here, we confirm that we can calculate the distribu¬ 
tion Vi{4>) of jammed packing fractions (for isotropic 
compression) using data from protocol 2. In Fig. [SJ we 
show that from protocol 1 agrees with Vi{(j)) = 

—MQ^dMi{(j))/d4> with 


Ml {(j)) = Mo exp 




(Al) 


of the simplest deformation methods is the application 
of simple shear strain 7 using Lees-Edwards boundary 
conditions. Thus, Lees-Edwards simple shear is strain- 
controlled. During the applied simple shear strain, one 
can measure the resulting stress of the system. 
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and the product given by the measured value of 


iM) - 


.F((^)5i(<^)< (</>)■ 


(A2) 


We find that (3 ~ 1.75 in 3D, whereas (3 ~ 2.0 in 2D. 


FIG. 6 : (a) Packing fraction 0 and (b) ratio of the stress 
anisotropy r to the pressure P versus shear strain 7 for all 
isostatic jammed N = 6 bidisperse disk packings. The solid 
black line in (a) obeys (f> = — 70 )^ -I- <j>o with A = 0.776, 

(fiQ = 0.665, and 70 = 0.35. Filled circles (downward trian¬ 
gles) indicate packings with positive (negative) local slope of 
0 versus 7 . 


Appendix B: Stress anisotropy in frictionless 
packings 

In this Appendix, we show the typical structure of the 
geometrical families for frictionless packings and measure 
the ratio of the stress anisotropy to the pressure for these 
packings. When deforming granular packings, the defor¬ 
mation method can be either strain- or stress-controlled. 
In strain-controlled deformations, a strain is applied to 
the system and the resulting stress is measured. In con¬ 
trast, in stress-controlled deformations, a stress is ap¬ 
plied to the system, and the resulting strain is measured. 
In simulations with periodic boundary conditions, one 


To measure the stress in 2D, we define the 2x2 stress 
tensor 

Ea 5 = I (Bl) 

i>j 

where fij\ is the A-component of the pairwise repul¬ 
sive force fij on particle i from particle j, Vijs is the 
d-component of the center-to-center distance vector rij 
between particles i and j, A = x,y, and S = x,y. In 
Fig. | 6 ] (b), we show the ratio of the stress anisotropy 
T = |Ei — S 2 I /2 to the pressure P = (Ei -|- S2)/2, where 
El and E 2 are the two eigenvalues of the stress tensor, 
as a function of the shear strain 7 . In many cases the 
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normalized stress anisotropy t/P follows nearly linear 
segments along the geometrical families, which appear 
as parabolas when 0 for each jammed packing is plot¬ 
ted versus 7 (Fig.|n](a)). However, for other geometrical 
families, r/P appears quadratic in 7 . For both cases, 
t/P decreases when —\d4)/d'^\ < 0 and increases when 
d(j)/d'y > 0 along each geometrical family. 


Appendix C: Comparison of pure and simple shear 

In this Appendix, we compare results for simple shear 
and pure shear protocols. All of the results for proto¬ 
col 2 presented in the main text were obtained using 
simple shear strain using Lees-Edwards boundary con¬ 
ditions. We have also studied strain-induced jamming 
using pure shear, where the separation between one pair 
of opposing edges of the simulation box is increased by 
1-I -7 and the separation between the other pair of oppos¬ 
ing edges is decreased by 1/(1 -I- 7 ). This deformation is 
the simplest example of a variable-shape simulation cell 
method that conserves volume. 



</> 


FIG. 7: Average shear strain 7j(</>) required to jam an ini¬ 
tially unjammed configuration at (f) using simple shear (light 
circles) and pure shear strain (dark squares) for N = 32 (blue 
symbols), 128 (green symbols), and 512 (red symbols). 


We show in Fig.[7]that the strain jj{(j)) required to jam 
an initially unjammed configuration at packing fraction 
(j) behaves qualitatively the same for packings generated 
via simple and pure shear strain. We see that the results 
for 7 j for simple and pure shear strain begin to devi¬ 
ate at low (j), but the deviation decreases with increasing 
system size. In the studies of pure shear, we stopped the 
simulations when the size of the cell in the thin direction 
caused interactions between a disk in the main cell and 
one of its own periodic images. If the system had not yet 
jammed, we did not include this trial in the measurement 


of 7 j. Thus, in the pure shear simulations, our results at 
low (j) and small system sizes were biased toward small 
strains. This effect vanishes in the large-system limit. 
Appendix D: Theoretical model 


In this Appendix, we elaborate some of the key aspects 
of the theoretical model described in the main text. We 
develop the theoretical model using an analogy with ab¬ 
sorption problems to calculate the probability to obtain 
isostatic jammed packings as a function of the path that 
the system traverses in configuration space. In princi¬ 
ple, the number density of jammed packings -F, depends 
on the 2N coordinates of the disks, the packing fraction 
(/, and boundary deformation 7 , but not the packing- 
generation protocol. After integrating over the 2N coor¬ 
dinates of the disks, is a function of and 7 . However, 
for the theoretical description, we assume that the num¬ 
ber density F'((/)) is only a function of packing fraction 
since P becomes independent of 7 in the large system 
limit as shown in Fig. 2 (c) in the main text. 

We imagine that a one-dimensional trajectory £{4), 7 ) 
through configuration space will encounter the basin 
of attraction of a jammed packing with a probability 
P{4>)S{4))dC during a step of size dC in configuration 
space, where S{(j)) is the 2A—1-dimensional cross-section 
of the basin of attraction of a jammed packing perpendic¬ 
ular to dC and dC = ^l{(j))d(j) and £ 2 (</')d 7 for protocols 
1 and 2 , respectively. 

For the results presented in this manuscript, our calcu¬ 
lations do not require complete enumeration of jammed 
packings and independent measurements of P{4>), S{(l>), 
and An advantage of our work is that we showed 

that one can obtain the product PS£ without complete 
enumeration by measuring the decrease in the number of 
unjammed configurations during shear. We showed (for 
fast quenching protocols) that the product PS£ depends 
only on the packing fraction (j), and not on the shear 
strain 7 . This result implies that we can use the shear 
protocol to predict the distribution of jammed packings 
obtained from the isotropic compression protocol. 

To obtain £ 1 ^ 2 , we measured the cumulative distance 
traveled by the system in configuration space after taking 
a step in packing fraction dcj) (Protocol 1) or a step in 
shear strain d'j (Protocol 2) and minimizing the total 
potential energy: 






N 

E 

2=1 


\Sri\ 


(Dl) 


where Sri is the change in position of particle i follow¬ 
ing the compression or shear step and subsequent energy 
minimization. We averaged ^ 1 , 2 ( 0 ) over at least 100 in¬ 
dependent trajectories. 
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